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The accelerating expansion of the Universe at recent epochs has called into question the validity of 
general relativity on cosmological scales. One probe of gravity is a comparison of expansion history 
of the Universe with the history of structure growth via gravitational instability: general relativity 
0^ ■ predicts a specific relation between these two observables. Here we show that the mean pairwise 

streaming velocity of galaxy clusters provides a useful method of constraining this relation. Galaxy 
cluster velocities can be measured via the kinetic Sunyaev-Zeldovich distortion of the cosmic mi- 
crowave background radiation; future surveys can provide large enough catalogs of cluster velocities 
to discriminate between general relativity and other proposed gravitational theories. 
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The most perplexing observation in physics today is the accelerating expansion of the Universe (for a review, see 
While such an acceleration can be brought about by a constant energy density of the vacuum, the associated 
energy scale is a small fraction of an electron volt. This energy scale is not, as far as we know, a natural fundamental 
u: scale in physics, and skepticism is warranted about new fundamental physics at room-temperature energy scales which 
^ r~| only manifests itself in cosmological phenomena. 

The standard hot big bang model of cosmology assumes that the Universe is statistically homogeneous and isotropic, 
and that its dynamics are determined by general relativity. The Einstein Equation describing the evolution of the 

■ metric then reduces to the Friedmann Equation, which can be written in the form 
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where p(t) and p(t) are the mean energy density and pressure of the Universe at a time t, and a(t) is the scale factor, 
giving the ratio of the separation between two objects in the cosmic rest frame at time t to their separation today. 
The scale factor is the function which describes the expansion history of the Universe, and the Friedmann Equation 
, is its dynamical equation. It is clear from Eq. (p} that, if general relativity is correct, we must have w = p/p < — 1/3 
for the expansion of the Universe to be speeding up at some particular epoch. Hypothetical stress-energy components 

|- | , obeying this relation have been termed "dark energy." Current observations show that w today is near — 1. 

The only logical alternative to dark energy which can explain the observational data is a modification of general 
Q\ , relativity, so that the dynamics of a(t) are determined by an equation different from Eq. ([1]). A variety of attempts 
have been made so far in this direction (see, e.g., @, S although modifying general relativity on cosmological 
J> , scales while still preserving successes on solar system scales and also matching available cosmological data on structure 
formation is challenging. Modified gravitation theories also tend to be more difficult to solve than general relativity, 
so detailed cosmological predictions for a given theory are often lacking. 



How can we distinguish between general relativity plus dark energy and modified gravity in a model-independent 
way? A number of papers have pointed out that in general relativity, a specific relationship exists between two basic 
gravitational phenomena in cosmology: the expansion history of the Universe and the growth of cosmic structure 
[a, Mi Hill- For scales well inside the cosmological horizon, the linear-theory growth factor D(a) is determined by the 
differential equation (see, e.g., [§]) 
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where the Hubble parameter H(a) = (1 / a) (da / dt) , Hq is the present value of H, and O m (a) = 8irGp m (a)/(3H(a) 2 ) is 
the ratio of the matter density p m to the critical density. (This equation assumes that any energy density components 
besides matter and radiation have negligible density variations.) The solution to this equation can be described to a 
very good approximation by 
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for a wide range of realistic models [5[ . The point behind this useful parameterization is the separation of the effect 
of expansion history, encapsulated in the function f2 m (a), from the growth rate, conveniently described by the single 
exponent 7. For standard cosmological models with dark energy, 7 ~ 0.55 + 0.05[1 + w(a = 0.5)] (reducing to the 
familiar 7 = 0.6 for w = 0). Sophisticated and general parameterizations of the evolution of the scale factor and 
structure formation in theories different than general relativity have been constructed 1^, but Eq. ([3]) provides 



a simple, single-parameter relation valid for general relativity that can be observationally tested in principle. While 
a given modified gravity theory is not guaranteed to have a substantially different value for 7 than general relativity 
with dark energy (see [12| for an example), generically this will be true. Linder and Kahn, for example, calculates 
that DGP gravity [l3( gives 7 ~ 0.68 [14| and give 7 for various scalar-tensor theories. 

Testing Eq. (|3|) observationally is not an easy task, however. It requires an observable which measures the growth 
rate of structure with good precision over a large range of redshifts. Directly measuring galaxy clustering as a function 
of redshift can in principle give the linear growth factor if the bias factor between galaxy clustering and the underlying 
mass distribution is understood well enough and if surveys complete to high enough redshift are available. Linder has 
advocated using redshift-space distortions, which are a measure of the internal velocity dispersions of bound objects 
like galaxy clusters [H| : up to a bias factor relating the galaxy distribution to the underlying mass distribution, the 
redshift space distortion is proportional to the left side of Eq. ([3]). This technique is promising, but relies on dynamics 
in the nonlinear regime, and requires spectroscopic redshifts of many galaxies. Both of these optical observation 
methods become increasingly difficult at high redshift. 

We advocate a different approach to testing Eq. ([3J: velocities of galaxy clusters obtained from the kinematic 
Sunyaev-Zeldovich (kSZ) effect. The microwave background radiation has its blackbody temperature shifted as it 
passes through a galaxy cluster, with the temperature shift being proportional to the line-of-sight velocity of the 
cluster and to its total optical depth for Compton scattering of the microwave radiation [l6| . For typical masses and 
velocities of large galaxy clusters, the temperature shift will be on the order of a few micro-Kelvin, on an angular 
scale of around one arcminute. The Sunyaev-Zeldovich effect is a powerful probe of cosmology because it is essentially 
independent of cluster distance. Also, the kSZ effect directly measures velocities with respect to the cosmic rest frame, 
unlike redshift-based velocity measurements which generally must contend with cosmological redshifts that are much 
larger than the redshift due to velocities. The kSZ effect, a temperature shift on the order of one part in a million 



of the mean background temperature, has not yet been detected (see 17|,[l8j for upper limits), but a new generation 



of experiments 1J3, 2^, 21 , 23, 23| are now making measurements at these angular scales and have the potential to 
detect the kSZ effect in clusters. 

Here we simply assume that future Sunyaev-Zeldovich surveys will result in a galaxy cluster velocity catalog, with 
each cluster having its sky position and redshift known exactly and its line-of-sight velocity determined with some 
characteristic error. From such a catalog, a number of different velocity statistics can be formed which are useful 
for cosmology. We have previously demonstrated the utility of cluster velocity statistics for constraining properties 
of dark energy [3, HE]- Consider in particular the mean pairwise cluster relative velocity 0, |27| . which can be 
estimated from only observed line-of-sight velocity components [25| . Using a pair conservation equation, an analytic 
approximation for the mean relative velocity for two clusters separated by a comoving distance r and at an average 
scale factor a is |28j 
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where £(r, a) is the cluster two-point correlation function and £(r, a) is the correlation function averaged over a sphere 
of radius r. Both of these correlation functions can be computed from the matter power spectrum and a bias giving the 
average number of clusters which form in a given overdensity. We assume that this bias is given by the standard LCDM 



cosmology 25[, which may not be precisely valid for alternate theories of gravity. However, numerical simulations 
[29l | suggest that deviations in the large-scale bias for alternate theories of gravity is only a few percent, so we expect 
this assumption to have little effect on our results. 

Notice that the amplitude of Eq. (j4|) is proportional to din D /din a, given by Eq. ((3j. So clearly this statistic can 
be used to measure the structure growth index 7. To quantify this assertion, we have computed the constraints on a 
5-parameter standard ACDM cosmological model from a cluster velocity catalog with a given number of clusters and 
a given mean velocity error, combined with priors on each parameter expected from the Planck Satellite's upcoming 
measurement of the primary microwave b ackg round temperature fluctuations [30| , and a measurement of Hq — 72 ± 8 
km/s/Mpc from the Hubble Key Project [3l| . We perform a standard Fisher Matrix estimate [12] of the constrained 
region in the multi-dimensional parameter space consisting of the amplitude of density fluctuations as , the power-law 
index of the primordial density perturbations n, the Hubble parameter today Hq, the present matter density fi m , 
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FIG. 1: The 1-sigma constraint on the growth index 7 in a seven-parameter spatially flat cosmological model, constrained by a 
galaxy cluster velocity survey of 4000 clusters chosen via the Sunyaev-Zeldovich effect, as a function of assumed velocity errors. 
Prior parameter constraints anticipated from the Planck satellite measurement of the primary microwave background power 
spectrum, plus the current Hubble Key Project constraint on the Hubble parameter, are included. 



and the growth index 7. For simplicity we assume the Universe is spatially flat, as indicated by current observations; 
the fiducial model for the scale factor evolution is standard ACDM. The resulting constraint on the growth index 7 
for a cluster velocity catalog with 4000 cluster velocities, marginalized over the other parameters, is shown in Figure 
1. The horizontal axis gives the standard error for measuring each cluster line-of-sight velocity and the vertical axis 
gives 17(7), the 1-er standard error on the resulting measurement of 7. 

For the current best-fit cosmological model, 4000 clusters corresponds to measuring a velocity for all clusters with 
masses larger than 2 x 10 14 M© in around 400 square degrees of sky. This mass is near the anticipated cluster detection 
threshold for current Sunyaev-Zeldovich experiments, although they are still a ways from having the sensitivity to 
measure the smaller kinematic SZ signal. If we measure the cluster velocities with an error of 400 km/s via their 
kinematic SZ distortion, this will provide a measurement of 7 to 0.08, a level which is interesting for discriminating 
various modified gravity theories from dark energy. 

This simple calculation likely gives a conservative estimate of 17(7), for a number of reasons. The error is statistics 
dominated, so it can be decreased by measuring velocities for clusters in a larger sky region, with 17(7) scaling like 
the inverse square root of the sky area. Measuring velocities more precisely can also somewhat increase the precision 
in measuring 7, although internal motions of cluster gas provide an astrophysical limit of around 100 km/s to how 



well cluster velocities can be measured with the kSZ effect 33j. The prior on cosmological parameters used here 
also does not account for correlations between various parameters constrained by the microwave background power 
spectrum; a more detailed parameter space investigation will likely result in smaller errors on 7. Including additional 
velocity statistics or other measures of structure formation also may decrease the error on 7, provided correlations 
between the various statistics are correctly accounted for. Finally, including 7 as an extra parameter generally does 
not substantially degrade the simultaneous constraints on other cosmological parameters [34[. 

The current experiments ACT and SPT are mapping large portions of the sky at arcminute resolution in multiple 
microwave frequency bands. Such measurements will constrain certain linear combinations of gas temperature, line- 
of-sight gas mass, and line-of-sight gas velocity, depending on frequencies and noise level 35(. The addition of a gas 
temperature determination from X-ray measurements often greatly improves the precision of cluster velocity deter- 
minations. While no cluster peculiar velocities have yet been measured, velocity catalogs for hundreds or thousands 
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of clusters are clearly within reach as the noise level of microwave maps decreases. 

The kinematic SZ effect for galaxy clusters provides a unique window into the growth of structure in the Universe. 
Like its thermal SZ counterpart, it is essentially independent of cluster distance, so it can probe structure growth over 
all epochs and over huge volumes. But in contrast to thermal SZ cluster detection to measure the evolution of cluster 
number density, cluster velocities derived from the kinematic SZ signal depend only weakly on cluster mass [36} (since 
the gravitational field causing cluster peculiar velocities provides the same acceleration to all masses), sidestepping 
systematic uncertainties related to the connection of cluster mass to observed SZ signal. The experimental challenge 
is daunting: detection of the tiny blackbody kinematic SZ distortion at arcminute resolution in multiple frequency 
bands, and disentangling this signal from other larger contributions including the thermal SZ distortion, dust emission, 
and sub-millimeter galaxy emission (e.g., §37} for a recent measurement). But progress has been rapid, and the payoff 
is one of the few reliable methods available for probing the fundamental properties of gravitation on cosmological 
scales, perhaps shedding light on the accelerating expansion of the Universe. 
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